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ABSTRACT 

We present deep HST/STIS coronagraphic images of the f3 Pic debris disk 
obtained at two epochs separated by 15 years. The new images and the re¬ 
reduction of the 1997 data provide the most sensitive and detailed views of the 
disk at optical wavelengths as well as the yet smallest inner working angle optical 
coronagraphic image of the disk. Our observations characterize the large-scale 
and inner-disk asymmetries and we identify multiple breaks in the disk radial 
surface brightness profile. We study in detail the radial and vertical disk structure 
and show that the disk is warped. We explore the disk at the location of the 
f3 Pic b super-jupiter and find that the disk surface brightness slope is continuous 
between O'/5 and 2'/0, arguing for no change at the separations where f3 Pic b 
orbits. 

The two epoch images constrain the disk’s surface brightness evolution on 
orbital and radiation pressure blow-out timescales. We place an upper limit of 
3% on the disk surface brightness change between 3-5”, including the locations 
of the disk warp, and the CO and dust clumps. 

We discuss the new observations in the context of high-resolution multi¬ 
wavelength images and divide the disk asymmetries in two groups: axisymmetric 
and non-axisymmetric. The axisymmetric structures (warp, large-scale butterfly, 
etc.) are consistent with disk structure models that include interactions of a plan- 
etesimal belt and a non-coplanar giant planet. The non-axisymmetric features, 
however, require a different explanation. 

Subject headings: planetary systems — planetary systems: formation — plane¬ 
tary systems: protoplanetary disks — stars: individual (Beta Pictoris) — minor 
planets, asteroids — techniques: high angular resolution 


Introduction 


Debris disks represent the late stage of planet formation when the primordial disk has 
dissipated and the accretion of giant planets has concluded, but rocky planets may still 
continue to accrete from a disk of icy/rocky minor bodi es. C ollisional cascades of the small 


bodies generate copious amounts of fine dust (e.g., Gaspar et al. 2012), often detectable 


m 


thermal emission and sometimes in scattered light. With typical ages between ten to a few 
hundred million years, bright debris disks provide insights into the properties and evolu¬ 
tion of young planetary systems. The distribution of fine dust in these systems can reveal 
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the s t ructure of the unde r lying planetary systems (e.g. iStark & Kuchnerl 120081: lApai et al. 
20081: iRoberge et al.ll2009l: ISu et al.ll2013l) and allows statistical st udies of the planetary sys¬ 
tem properties, dynamics, a nd collisional history (for reviews see iMever et al.l 120071: IWvatt 
20081: Matthews et al.ll2014p) . In addition, debris disk s likely pose a critical limitation on fu¬ 
ture missions to directly image terrestrial exoplanets (iGuvon et al.ll2006l : iMillan-Gabet et al. 

20 Ilk 


However, debris disks around stars known to host exoplanets remain very rare, making it 
difficult to establish a quantitative link between disk structures and planets. f3 Pic is a unique 
system hosting the best-studied debris disk and in it_a directly imaged super-jupiter. This 
system is also nearby (19.44±0.05 pc, van Leeuwen 20071) . allowing high physical resolution 
with current instrumentation. 


The [3 Pic disk is massive and exceptionally bright {F^^JF+ = ^2.5 x 10 , lLagrange et al 


2000 1). with a diameter exceeding 400 au (e.g. Smit h k Tcrrilc 19841 ). The disk has been 
imaged with multiple instruments, but the most specta cular and most detailed images have 


been obtained with HST/STIS and HST/ACS/HRC (iHeap et al.l l2000l: iGolimowski et al. 


20 06). Perhaps the most perplexing structure in (3 Pic is a warp or secondary disk inclined 
~ 5° to the main disk. The STIS coro nagrap hic images have been interpreted as a warped 


disk (Heap et al.l 2000), while Golimowski et al. (120061) have used deconvolved HST/ACS 
images to argue for the presence of a primary and a smaller, fainter, and inclined secondary 
disk. In the deconvolved ACS images the two wings of the main disk show prominent NW-SE 
asymmetry and at least three breaks between 40 and 250 au in its radial surface brightness 
distribution, indicative of changes in the disk structure or in dust grain properties. These 
breaks in surface bright ness distribution may emerge from dust belts and/or planetesimals 


belts (e.g. Wilner e_t ah 20111). The NW and SW wings of the disk also have different optical 


colors; in addition, while the NW wing is linear, the SW wing is bowed. 

I 11 the HST/ACS images the two wings of the secondary disk follow similar radial bright- 
ness p rofile s, su ggesting identical radial density distributions and dust grain populations 
(Golimowski et al. 20061) . Surprisingly, however, these two linear wings are not collinear. In 
summary, the scattered-light f3 Pic disk has a complex structure with evidence for different 
grain populations and radial density profiles, as well as breaks, discontinuities or asymme¬ 
tries in the central region of the disk. Planetary perturbations have been inv oked to ex plain 
the large-scale asymmetries in the main disk or the secondary disk/warp (e.g. iMouillet et al. 
1997 : Augereau et ah 2001 : Golimowski et al. 2006 ). 


Recently, lLagrange et al.l (120091) and lLagrange et al.l (120101) announced the discovery 
and confirmation of a giant planet orbiting j3 Pictoris. The planet has an estimated mass 
of 8-10 Mj and an orbital semi-major axis of about 8 au (Chauvin et al. 2012). Since the 
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discovery of the planet its motion has been closely monitored and its orbit has been gradually 


refined (e.g. C 

)uanz et al. 

2010; 

Bonnefov et al. 

2011; 

Chauvin et al. 

2012; 

Males et al. 

2014; 

Nielsen et al. 

2014; 

Macintosh et al. 

2014 

). The planet’s orbit was found to be slightly 


inclined 1°) with respec t to the main disk an d to have a relatively low eccentricity 


(e<0.1, lLagrange et al.l 12012c Macintosh et al.ll2014Jf . 


Recent observations with Herschel (jVandenbussche et al.ll2010l) and SMA (IWilner et ah 


2 011 ) sampled the spectral energy distribution of /3 Pic at far-infrared and millimeter wave¬ 
lengths and provided resolved or marginally resolved images of the disk. Most interestingly, 
Wilner et al. (2 0111) identified two millimeter clumps at radii ~375 and argued that these 
indicate a population of planetesimals in a ring at r = 94 ±8 au - inclined with re spect to t he 


midplane - whose collisions serve as a source of dust grains. A few months ago Dent et al. 


(2014) published ALMA high-resolution (12 au) and high sensitivity 870 /ini continuum and 
12 CO 3-2 transition observations of the /3 Pic disk. The continuum images revealed that the 
southwestern disk is brighter than the northeastern one and shows a peak at 60 au (3708) 
with a tail extending beyond 470. The CO distribution also shows excess emission in the 
SW side, but it peaks at a larger separation (85 au or 474). Unlike the continuum emission, 
the CO clump detected by ALMA peaks above the disk midplane (~5 au). The velocity- 
based de-projection of the CO data argues for, but does not prove, a broad CO clump at 
the Earth-facing SW quarter of the disk at with an orbital radius of 85 au and a CO tail 
extending toward the NE side of the disk (iDent et al.l 2014). The dust mass is estimated 
to be 4.7 ± 0.5 x 10 23 kg (6.4 M Moon , while the CO mass is approximately 1.7 x 10 2 ° kg 
(0.0023 M Moon ). The presence of such large amounts of CO gas is surprising, given its short 
dissociation lifetime (~120 yrs) and poses a problem similar to that seen in an another old 
debris disk (HD 21997, Kospal et al. 2013 ). 


D ent et aL| (2014) identify two possible scenarios to explain the morphology of the CO 


and continuum observations: In the Erst, the clumps emerge from collisions of planetesimals 
trapped in 2:1 and 3:2 mean motion resonance with an outward migrating >10 M Earth 
planet (W yattl 120031) . In the second, a single, massive, and recent (~0.5 Myr) collision with 


an approximately Mars-sized parent body injects dust and CO gas into the system. 


In a 2014 VLT adaptive optics study Milli et, al. (2014) imaged the disk between 074 
and 378 in L'. They report an overall bowed disk structure, particularly prominent in the 
inner disk (<1” from the star). Based on an anisotropic scattering model they explain the 
key features in the disk surface brightness distribution (bowed structure, surface brightness 
profile, and the warp) by a two-component disk: an outer (main) disk extending outward 
from an inner radius ( r warp ), and an inner warped disk inclined at 4° with respect to the 
main disk. In their model, the entire disk is inclined toward the observer at an inclination 
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angle 86°. 

In this study we explored two questions important to understanding the /3 Pictoris 
system: i) What is the structure of the inner disk?, and ii) Are there planetesimal groups 
in orbits resonant or co-moving with f3 Pic b? We addressed these questions by obtaining 
high-quality coronagraphic images of the system using the Hubble Space Telescope. The new 
optical disk images have the smallest yet inner working angles (~ (y/35), allowing detailed 
studies of the inner disk in the absence of image artifacts at small stellocentric distances 
previously probed. We also re-reduced data from 1997, obtained in a similar mode but 
with a larger inner working angle, (now) with PSF-template subtracted coronagraphy. By 
comparing the disk structure in the two epochs we searched for temporal evolution in the 
inner disk over timescales comparable to the orbital period of the planet and similar to 
radiation pressure timescales. 

The paper is organized as follows. We first review the observations, data reduction 
and basic data analysis (Section [2]). Section l3Tl describes the disk’s orientation and overall 
shape, Section 13.21 discusses the surface brightness profiles along various directions and the 
large-scale asymmetry of the disk, Section 0] explores the vertical structure of the disk and 
the warp through various methods, Section [5] compares the orbit of the planet to the inner 
disk structure and explores temporal evolution of the scattered light disk. Section [6] provides 
a multi-wavelength view of the disk. Section [T] discusses our observational results in the 
context of disk-planet interaction models. Finally, Section |S] summarizes the key findings of 
our work. 


Observations and Data Reduction 


In this paper we pr ese nt a re-reduction of a 1997 HST/STIS (Space Telescope Imag¬ 
ing Spectrograph, Woodgate et ah 19981) archival dataset and also introduce a new set of 
HST/STIS coronagraphic observations we obtained in 2012. STIS coronagraphy utilizes an 
image plane mask including two orthogonal wedge occulters that is inserted in the STIS fo¬ 
cal plane. STIS cannot use filters in its coronagraphic (50CORON) mode; thus, the images 
described below are unfiltered data whose wavelength coverage is only set by the spectral 
response of the STIS CCD detector (~200-l,050 nm with A pi „ ot =575.2 nm and FWHM=433 
nm and a spectrally flat source). Based on the STIS Exposure Time Calculator 22.1.1 and 
assuming an input stellar spectrum of A5V (corresponding to f3 Pic ) we End that the effec¬ 
tive wavelength (the wavelength at which incoming photons generate the most electrons) is 
571.3 nm. The STIS CCD detector is a 1024x1024 pixel UV-enhanced chip with an image 
pixel scale of 0 , . , 05077 per pixel and at these wavelengths a resolution element is ~60 mas. 
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2.1. The 1997 STIS Observations 


The 1997 observations and their analysis are described in detail in Heap et ah (2000) 
and we only give a brief summary of the data we have used here. Table [Tj summarizes the 
observations.The first epoch observations were carried out in program GO-7125 (PI: Heap) on 
Sep 16-17, 1997 with STIS occulting wedge B. (3 Pic was observed in three non-contiguous 
orbits (program visits #4-6, archived raw dataset IDs: 04204, 05, 06). Prior data from 
visits 1-3 using wedge A failed or were degraded and ill-suited for an investigation. In the 
first wedge B orbit the disk was placed halfway between the telescope’s diffraction spikes 
and (nearly) orthogonal to the occulting wedge, while the second and third orbits were 
executed with spacecraft off-rolled by —12° and +14° relative to the first orbit. During the 
observations /3 Pic was occulted by WedgeB2.0 and WedgeBl.O. At each location a series of 
images were taken with 5 s and 3 s integration times, respectively. Here we only re-reduced 
the WedgeBl.O images, for lack of a suitable Wedge B2.0 reference PSF. 


2.2. The 2012 STIS Observations 


Our new STIS observations were carried out on March 6, 2012 and used a refined 
observing strategy that resulted in higher signal-to-noise, we aker PSF residuals, and even 
smaller inner working angle. Our observing strategy followed Schneider et al. (2014), using 
multiple roll angles, coronographic wedge positions, and contemporaneous observations of a 
color- and brightness-matched reference star to minimize PSF residuals. 


The data were acquired in program GO-12551 (PI: Apai). Table Q] gives a summary 
of the observations. In this program using three contiguous orbits we targeted [3 Pic in 
the first and third orbits and the PSF reference star (a Pic) in the second orbit. The 
PSF star (a Pic, V=3.3, B—V= 0.18) was selected to closely match the celestial position, 
apparent magnitude, and color of (3 Pic (V=3.86, B—V=0.17). We executed spacecraft rolls 
between the three orbits to recover (in combination) azimuthal sectors in the celestial frame 
otherwise unimaged due to the imposition of the coronagraphic wedge or degraded by the 
telescope diffraction spikes, and to further reduce the rotationally non-invarient components 
of the PSF-subtraction residuals. In each of the three orbits we used a combination of short, 
medium, and long integrations (13.2 s, 48 s, 240 s) at 0'.'6 and l'.'O wedge width locations; 
while the 0'.'6 location was observed in WedgeA, the l'.'O wedge location was observed both 
at the A- and B coronagraphic wedges of STIS. The inner working angle (IWA) equals the 
half-width of the occulting wedge at the stellar position: WEDGE(A or B)0.6 = 0'.'3 IWA 
and WEDGEA—B1.0=0'.'5 IWA. Due to the very high contrast in the image between the 
star, inner disk, and the outer disk no single image has suitable dynamic range to capture 
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the scattered light surface brightness levels at all radii; therefore, we used short integrations 
to probe the stellar PSF and the inner disk, and longer integration times to probe the 
outer disk (while saturating the inner disk). The combination of the non-saturated pixels 
in images with different integration times allowed us to significantly increase the effective 
dynamic range of our final images and to probe the disk with high signal-to-noise over a 
surface brightness range of more than 10 5 . For more details on the image combination we 
refer the readers to (1 Schneider et ah 2014). 


2.3. Coronagraphic Data Reduction 


Our data reduction follows that of Schne ider et ah (120091 1. ffere we provide a brief 
summary of the key steps and differences from that dataset and note that the GO 7125 and 
GJ12551 images were reduced separately. 

First, the RAW hies were reprocessed with the STIS basic calibration package calstis 
using the most up-to-date dark and bias reference hies closely contemporaneous with the ob¬ 
servations. Each instrumentally calibrated image was then manually inspected for anomalies, 
but none were found. Within each visit groups of exposures with the same setup (target, 
orientation, pointing, aperture, and exposure time) were median combined. Next, for each 
median-combined image we used small sub-arrays located far from the target to estimate 
the combination of the sky and instrumental background, which were then subtracted from 
the images. By dividing each background-subtracted median combined image with the ex¬ 
posure of the individual images we calculated the count rates per pixel. We next located the 
positions of the targe t (oc cul ted star) by identifying the intersection of two lines ht to the 


diffraction spikes (see [Sc h ne ider et al. 2014 for details of the process). 


The above steps have been repeated on every frame taken on f3 Pic and the PSF refer¬ 
ence star. In the next steps from each (3 Pic image we subtracted a matching PSF template. 
Because in the GO 7125 no PSF template star observations were taken, from these images we 
subtracted the WedgeB PSF template derived from the 2012 epoch GO 12551 observations. 
The precise alignments of the target and PSF images are critical and we used a two-step 
alignment procedure. First, the target and template images have been co-aligned using frac¬ 
tional pixel o ffsets to the initial centers measured from the diffraction spike-fitting procedure 


([Schne i de r et al. 2009J, §4.2). Then we refined the alignments between the target and PSF 


images by iteratively minimizing the diffraction spike residuals in the PSF-subtracted im¬ 
ages where the disk flux did not dominate by varying three free parameters (A X and A Y, 
and intensity scaling). We iterated all three parameters until convergence was reached in 
minimizing the subtraction residuals. This solution was also verified by visual inspection. 














Table 1: Summary of the STIS observations used in this work. 


Prog. 

Visit # 

Target 

Date 

Int. Times [c 

5 

ORIENT 61 

Aperture 

7125 

04 

f3 Pic 

1997-09-17 

16x3 


210.07 

WEDGEB1.0 

7125 

05 

/3 Pic 

1997-09-17 

16x3 


224.07 

WEDGEB1.0 

7125 

06 

/3 Pic 

1997-09-17 

16x3 


198.07 

WEDGEB1.0 

12551 

01 

f3 Pic 

2012-03-6 

11x1.2 


239.15 

WEDGEA0.6 

12551 

01 

13 Pic 

2012-03-6 

16x3.0, 4x60.0 

239.15 

WEDGEA1.0 

12551 

01 

/3 Pic 

2012-03-6 

4x60.0, 16x3.0 

239.15 

WEDGEB1.0 

12551 

02 

a Pic 

2012-03-6 

11x0.7 


245.96 

WEDGEA0.6 

12551 

02 

a Pic 

2012-03-6 

4x36.0, 16x1.9 

245.96 

WEDGEA1.0 

12551 

02 

a Pic 

2012-03-6 

16x1.9, 4x36.0 

245.96 

WEDGEB1.0 

12551 

02 

a Pic 

2012-03-6 

11x0.7, 2x(3x0.7) fe 

245.96 

WEDGEB0.6 

12551 

03 

f3 Pic 

2012-03-6 

11x1.2 


272.62 

WEDGEA0.6 

12551 

03 

(3 Pic 

2012-03-6 

16x3.0, 4x60.0 

272.62 

WEDGEA1.0 

12551 

03 

/3 Pic 

2012-03-6 

4x60, 16x3.0 

272.62 

WEDGEB1.0 


“ORIENTAT: position angle of the image +Y axis measured eastward from celestial north. 
^Additional two sets of three exposures each with | pixel in Y offsets for calibration 

purposes. 
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Fig. 1.— Direct comparison of the two STIS data sets on the Pictoris disk. Left: Onr 
re-reduction of the STIS image taken in GO-7125 in 1997. Right: New STIS images taken in 
our program GO-12551 provide higher signal-to-noise and better PSF-subtraction, as well as 
smaller inner working angle. The two images are shown here in a logarithmic stretch, with 
instrumental brightness units of counts per second per pixel; North is up. 
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This two-step alignment procedure was repeated for every target-PSF image pair. 

After precise target and template alignment and PSF subtraction the images were ro¬ 
tated to a common ” north up” orientation about the occulted star; image edge clipping was 
avoided by embedding the images in a larger, zero-padded image. For each image we created 
a binary ’’bad” data mask that flagged all pixels that were saturated, corrupted by diffrac¬ 
tion spikes, covered or affected by the wedges, or otherwise significantly degraded. For the 
short Wedge0.6A/B images this mask also included pixels at large stellocentric angles where 
the signal-to-noise was low. Finally, the rotated, masked images were median combined to 
create the final analysis-quality images from the GO 7125 and GO 12551 data sets. 

Our final images reveal the complex debris disk around (3 Pic between 0'.'35 to about 
13”, providing an image with the highest quality yet and also the smallest working angle 
(Fig. IT]) at visible wavelengths. 

Figure [3] shows our estimated signal-to-noise maps of the re-reduced 1997 and the new 
2012 images. The maps were calculated by dividing the instrumental counts by the noise 
in the images, both averaged over 3x3 STIS pixels. The noise for each pixel position was 
estimated as the standard deviation of the counts in the given pixel in each of the frames 
that covered that pixel. Note, that the reliability of the noise estimate for any given pixel is 
sensitive to the number of valid pixels that covered that location. 

To demonstrate the improvement in the image quality with the GO 12551 data to the 
re-reduced GO-7125 data we note that the signal-to-noise integrated over 2x2 pixel area in 
the 1997 images is 20, while it is 40 in the 2012 images (in the disk midplane at separation 
35 pixels or 1'.'77); in the better-sampled closer in regions the improvement is even more 
significant: the 2012 images reach a signal-to-noise of 120, while the 1997 images have a 
signal-to-noise of 12 (at 21 pixels or l'/06) after our reprocessing. In Figure[5]we plot the two 
images with the main disk axis aligned horizontally (assuming a position angle P.A.=29°.l) 
and by applying a radial brightness scaling (xr LS ) that allows easy comparison of the inner 
and outer disk structures. Here the 1.8 exponent is not physically motivated but was chosen 
for visual illustration in collapsing the radial dynamic display range. 


3. Disk Orientation and Surface Brightness Profile 
3.1. Disk Position Angle 

We determined the main disk’s position angle by measuring the directions of vectors 
pointing from the star’s position to the disk isophotically-determined mid-plane at stcllocen- 
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Radially-weighted Direct Comparison of 1997 and 2012 Images 



-15 -10 -5 0 5 10 15 

Distance ["] 


Fig. 2.— Direct comparison of the two disk images obtained from the re-reduction of the 
1997 dataset and from the new 2012 STIS images. While the 1997 data has a slightly larger 
effective field of view (not shown here), importantly the 2012 data have smaller inner working 
angle and higher signal-to-noise. These images have been multiplied by an r 1 ' 8 function to 
enhance the fainter outer disk’s visibility while simultaneously showing the inner disk. The 
image has been rotated counter-clockwise by 60.9°; the NE side of the disk is on the left. 
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Signal-to-Noise Map for 2012 Image 
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Signal-to-Noise Map for 1997 Image 
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Fig. 3.— Signal-to-noise ratio maps for the combined 2012 and 1997 images. The 2012 
images provide higher signal-to-noise ratio across the image and, in particular, at small 
inner working angles close to the star. The images have been rotated counter-clockwise by 
60.9°; the NE side of the disk is on the left. 
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trie separations of 11.2” in both the northeastern and southwestern wings for both the 1997 
and the 2012 images. Because the disk is not completely straight we did not attempt an 
algorithmic determination of the disk angle but instead opted for the manual identification 
of the brightest pixel at the given radii and calculated the slope of the line connecting these 
points. With this procedure we found a disk position angle of 29.17° for the 1997 image 
and a disk position angle of 29.05° for the 2012 image, each with an estimated uncertainty 
of ±0.1°. This uncertainty corresponds to 1 pixel uncertainty in the position of the disk 
midplane over the full visible disk in our images, which is realistic given our manual fitting 
procedure. However, this uncertainty does not i nclud e the uncertainty in the spacecraft roll 
angle, which is typically less than 0.1° fjKalas et al. 120131 ). To estimate the uncertainties 
in the absolute spacecraft roll angle between the two epochs we searched for the minimum 
of the difference between the x-y aligned 1997 and 2012 images on a finely sampled grid of 
rotation (0.05°/grid point). We found that a 0.0° relative rotation produces the smallest 
difference, i.e. the HST astrometric reference frame is highly reproducible between the two 
images, even though different guide stars were used in the two epochs at different spacecraft 
orientation angles. This finding reinforces our conclusion that the uncertainty of our disk 
angle measurement is dominated by the uncertainty in the disk midplane determination , 
and not differences in the spacecraft’s astrometric reference frame. Therefore, we adopt the 
average of the two position angles as the disk position angle, i.e. we will use 29.1±0.1° 
(counter-clockwise from north). 

Recently, it was pointed out that the /3 Pic disk midplane is not straight, but slightly 
bowed (Milli ct ah 2014). We test evidence for the overall bow in our data by independently 
fitting the two sides of the disk. We used the analysis quality reduced images from the 2012 
observations and fit two lines starting from the position of the star to the brightest point in 
the disk at the edge of the field of view (~11” on each side, same points as used above for 
the disk position angle determination, see black lines in Fig. UJ). These fits yield about 0.23° 
difference between the position angle of the two sides. As illustrated in the inset of Fig. Efl 
the star’s position lies NE from the vector connecting the SW and NE disk midplane (white 
line), which is consistent with a curvature to the NW direction, i.e. the two disk wings 
slightly bowed to the SE direction. In other words - with respect to a straight line fitted on 
the disk - the inner disk appears slightly NW of the line and the two disk wings (in NE and 
SW directions) both are slightly bent to the other sid e (SE). The fact that the disk is slightly 
bowed may be due to the interactions with ISM (e.g. iGaspar et al.ll2008l : IBnenzli et al .11201 f)h 
or, more likely, due to the effect of forward-scattering grains in a not completely edge-on disk 
fiRodigas et al. 2014); given the direction of the curvature and the inner disk asymmetries, 
it is not likely that the overall bow is a result of radiation pressure-driven grains emerging 
from the asymmetric inner disk. 
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Fig. 4.— Our 2012 image and the illustration of the disk angle determination. The disk 
position angle was determined by fitting a line (white line) on the brightest points in the 
disk (disk spine) 11” on either sides of the disk. This line does not precisely intersect the 
star itself (see inset), demonstrating an overall bow in the disk. The black lines connect the 
star to the brightest points on either side of the disk at 11” separation. The same procedure 
was repeated for the 1997 image, which gave very similar results. In this figure north is up. 


Disk Geometry and Position Angles 



-10 -5 0 5 10 

Arcseconds 









15 


3.2. Surface Brightness Along the Disk Midplane 


To measure the disk surface brightness profile we first rotate the ’’north up”, PSF- 
subtracted image of the disk (by 60.9° CCW) to place the main axis of the disk on the 
image horizontal and then extracted a 4 pixel-wide rectangular stripe centered on the disk 
mid-plane. The image was converted to physical flux density units by multiplying by the 
STIS inverse sensitivity (4.0169 xlO _19 ergs/s/cm 2 /A per count/s for the 1997 image and 
4.1446 x 10~ 19 ergs/s/cm 2 /A per count/s for the 2012 image), a value that is provided by 
STScI’s and is based on the absolute flux calibration of the STIS instrument. Next, for each 
radius (at each pixel) we calculated the median of the vertical stripe, excluding pixels not 
sampled in the image. 


We applied this process to the northeastern and southwestern wings of the disk and for 
both epochs (1997 and 2012), resulting in four surface brightness profiles. Figure [5] shows the 
resulting surface brightness profiles for the 1997 and 2012 data.The 2012 STIS observations 
have closer inner working angle and probe the disk surface brightness all the way down to 
0735; however, we consider the data reliable (i.e. sampled by multiple pixels at a given 
radius) only for radii greater than 074. 

In Fig. Owe also include representative ±1 a uncertainties for two of the profiles as blue 
and red shaded envelopes. We estimated the uncertainties at each radial bin as the standard 
deviation within the 4-pixel wide stripe, which reflects the combination of photon and read¬ 
out noise (minor contributions), PSF subtraction residuals, and actual vertical differences 
in the surface brightness in the disk (major contribution). The fact that the actual scatter 
in our profiles is smaller than the envelope suggests that our uncertainty estimate is conser¬ 
vative. We also considered another component, HST’s photometric calibration uncertainty, 
but concluded that it is not significant for our analysis for the following reasons. The cali¬ 
bration uncertainty is less than 1% (see §, 16.2[i and would therefore not visibly change our 
current noise estimate. In addition, it would effect all measurements the same way (small 
relative increase in surface flux density), i.e. and would therefore not change our results on 
the surface brightness slopes (see below). 


In Fig. [5] we also plot in light gray the ALMA 1.3 mm dust continuum profile (integrated 
in the direction perpendicular to the disk midplane) from the measurements presented in 
Dent et ah (2014). The nun-sized grain population traced by the ALMA observations is 


thought to follow the distribution of the planetesimals in the disk, but it is somewhat com¬ 
plicated by the projection effect in the nearly edge-on viewing geometry. We point out that 
the regions R2-R6 - covering annuli between break points in the surface brightness slopes - 
appear to all coincide with inflection points in the dust continuum. This agreement suggests 
that the fine grains seen in scattered light provide also approximately trace the distribution 
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of planetesimals. 


3.3. Large-scale Radial Asymmetry 

We confirm a previously reported brightness difference between the NE and SW wings 
of the disk. The SW side of the disk is brighter in the inner 8'ZO and fainter at larger distances 
on the NE side (see Fig. [5]). We find that this brightness difference extends to the previously 
unexplored inner disk down to at least 075 or about 10 au projected radius. The brightness 
difference is greater at the smallest separations (~50%), initially decreases toward larger 
radii, and the two sides reach identical brightness at 275. Between 370 and 870 the SW side 
is again brighter, but its brightness drops faster beyond 770 than that of the NE side and 
beyond 870 the NE side is seen to be brighter. 

We note that the more precise characterization of the surface brightness profile is limited 
by the fact that the disk has considerable vertical structure (§ 2J) and thus projecting it to 
a one-dimensional distribution is an ill-defined problem: different vertical aperture widths 
will return slightly different results depending on what disk regions are included. 


3.4. Breaks in the Radial Surface Brightness Profile 

The radial surface brightness profiles show distinct disk regions with different slopes. 
To characterize the surface brightness profiles we fitted linear slopes to the logarithms of the 
projected separation and the surface brightness profiles (see Fig. [5]). We find that the slopes 
are very well constrained and very similar in the two epochs. We note that the uncertainty 
of the slopes is not dominated by noise in the images but by the precise definition of the 
aperture over which the slopes are calculated in the non-linear and extended disk. We 
estimate the typical uncertainty of the power-law coefficients to be ±0.04. 

Our analysis of the disk’s radial surface brightness reveals at least five distinct disk 
regions (see Table [2]) ranging from 074 to 1374, corresponding to 7.7 au to 260 au. We 
labeled these R2-R6 and also included region Rl, which is only probed in our 2012 image. 
We characterized the surface brightness slopes by fitting the brightness distributions in log- 
log space with linear functions, corresponding to S oc r a power-law radial distributions. The 
surface brightness slopes determined range from a mild slope in the inner disk (a ~ —1.7 to 
— 1.8 for radii smaller than 39 au, i.e. Regions Rl and R2), surrounded by an even milder 
slope (a ~ —1.4 for radii between 39 au and 68 au, Region R3). Beyond 70 au we find three 
disk regions with much steeper slopes: the first has a ~ —1.9 (between 72 au and 110 au, 


17 


Disk Radial Surface Brightness 



Distance farcsecj 

Fig. 5.— Comparison of the disk surface brightness profiles from two epochs and two sides of 
the disk. The disk shows an NE-SW asymmetry: the NE (red) side is slightly fainter in the 
inner 8 r /0, but brighter at larger separations than the SW (blue) side. Multiple breaks are 
present in the disk surface brightness: at 2'.'0 in the SW wing only and at b'/O in both wings. 
The profiles observed in the two epochs show near-perfect matches, except for a possible 
fainter emission in the newer (2012, red dashed line) NE measurements at radii between l'.'O 
and 1"7”. change The light gray curve shows the ALMA 1.3mm continuum profile (Dent 
et al. 2014) integrated perpendicular to the disk and presented here in a linear, normalized 
scale. Note, that several of the disk region boundaries identified in the scattered light image 
the correspond to inflection points in the planetesimal distribution traced by the ALMA 
continuum data. The shaded region marks the probable semi-major axis range of /3 Pic b. 
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Region R4), the following large region shows the steepest slope observed (a ~ —4.4 for radii 
between 130 au and 190 au, Region R5). The outermost disk region again shows a steep, 
but somewhat milder slope (a ~ —3.9 to —4.5 for radii between 190 au and 260 au, Region 
R6). 


The disk surface brightness profile has been studied in detail by Go l imowsk i et al . (l2006i) 
on the basis of their HST/ACS F606W coronagraphic observations. The boundaries of our 


regions R3-R6 were matched to those of the four regions identified by Goli mowsk i et ah 


( 20061 ) (regions 1-4 in their paper) allowing dir ect co m pari son (see Column 4 in Table |2]). 
We note here that we adopted the slopes Golimowski et ah (20061) derived from the non 
PSF-deconvolved ACS images; these sho w greater similarity to our ST1S images than the 
PSF-deconvolved ACS images of Golimowski et ah (2006) . We find that the overall radial 
disk structure observed in our STIS images matches closely that described by those authors. 
In spite of the good overall agreement the surface brightness slopes fitted to the two different 
datasets are slightly different, which we attribute to the different spectral coverage of the 
unfiltered STIS observations (~200-l,050 A) and the filtered F606W ACS images. 


Our 2012 images have smaller inner working angles than the previous ACS and STIS 
observations, allowing us to accurately measure the radial surface brightness of the inner 
disk. In the 2012 data set we can probe the disk structure down to 0/4. We find that the 
slope observed between O'/S and 2/0 in the SW disk continues inward in the NE disk without 
significant change down to radii 0/4. 


4. Vertical Disk Structure and Asymmetries 


first imaged (e.g. 


Smith & Terrile 

1984; 

Kalas & Jewitt 

1995; 

Heap et al. 2000; 

Golimowski et al. 


20061 ). Our high-quality and small working angle images allow detailed study of the inner 
disk structure. In the following we use different analysis techniques to probe the asymme¬ 
tries in the disk. First, we show and discuss a radially normalized disk image and model the 
vertical structure at each radius as a sum of two gaussian functions. Second, we will explore 
the rotationally asymmetric disk component. Third, we attempt to separate the main disk 
contribution from the warped component. 
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Table 2: Power-law fits to radial surface brightness profiles and comparison to the values 
measured by Golimowski et al. (2006) (see their Table 3). 


Region 

Range 

Range 


a 



[arcsec] 

[au] 

1997 STIS 

2012 STIS 

F606W ACS a 

Northeastern Disk 

Rl 

0.4-0.8 

7.8-16 

-1.17 ±0.43 

-1.80 ±0.20 

— 

R2 

0.8-2.0 

16-39 

-1.80 ±0.05 

-1.87 ±0.10 

— 

R3 

2.0-3.5 

39 - 68 

-1.36 ±0.06 

-1.26 ±0.07 

-1.34 ±0.02 

R4 

3.7-5.6 

72 - 110 

-2.05 ±0.06 

-1.99 ±0.07 

—1.91±0.02 

R5 

6.6-10.0 

130-190 

-4.23 ±0.07 

-4.19 ±0.08 

—4.19±0.02 

R6 

10.0-13.4 

190 - 260 

-4.25 ±0.57 

—3.92 ± 0.49 b 

—3.63±0.02 

Southwestern Disk 

Rl 

0.4-0.8 

7.8-16 

-1.91 ± 1.63 

-1.58 ±0.33 

— 

R2 

0.8-2.0 

16-39 

-2.00 ±0.11 

-1.86 ±0.11 

— 

R3 

2.0-3.5 

39-68 

-1.45 ±0.16 

-1.37 ±0.16 

—1.06±0.02 

R4 

3.7-5.6 

72-110 

-1.89 ±0.18 

-1.85 ±0.18 

—2.03±0.02 

R5 

6.6-10.0 

130-190 

-4.83 ±0.09 

-4.70 ± 0.09 

—4.76±0.02 

R6 

10.0-13.4 

190-260 

-5.00 ±0.52 

—4.36 ± 0.63 b 

—4.00±0.02 


“The HST/ACS uncertainties may not include the systematics fully. 
& The outer radius of Region R6 in the STIS 2012 images is ll'/O. 
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4.1. Radially Normalized Vertical Structure 


We now inspect the disk vertical structure by modeling it with analytical functions at 
each radius. To better visualize the vertical structure we normalize the disk brightness at 
each radius (1 pixel vertical slice) by the peak value at each corresponding radius (see panel 
with Radially Normalized Image in Fig. EH and a 4x vertically stretched version in Fig. [7]). 
This image reveals a prominent asymmetry in the disk: the disk is warped in its NE (left 
in Fig. [6]) wing between ~ 2'! and ~ 5”; a similar, but not identical warp is apparent in 
the SW wing between ~ 3"0 and S'.'S. The warp, in effect, introduces a counter-clockwise 
deformation in the disk, i.e. adding excess scattered light south of the NE wing and north 
of the SW wing. 


To quantify the vertical structure we fitted each vertical slice with the combination of 
two Gaussian functions. In Fig. [6] we fixed the center of one of the Gaussians to the disk 
mid plane, while in Fig. [7] we treated the centers of bot h G aussia ns as free parameters. The 
fits were optimized by the IDL function MPFIT (Markwardt 2009) which applies Levenberg- 
Marquardt least-squares optimization. We found that the sum of two Gaussians provided a 
very good fit for all radii, with the exception of the innermost 1”. The top panels in Figs. [6] 
and [7] show the vertical disk brightness profiles (blue), the best-fit analytical functions (red), 
and the residuals (thin black lines) for the four vertical slices shown in the radially normalized 
image (labeled as Cuts A-D). 


These best-fit analytical functions allow us to generate a simple ” disk model”, which we 
subtract from the actual radially normalized image to verify how well this simple approach 
can reproduce the disk structure. In the second panel from the bottom of Fig. 12 we show 
the model-subtracted radially normalized disk. This image demonstrates that our simple 
analytical fit captures well most of the disk structure down to a stellocentric separation 
of 1”. We note that the fit shown in Fig. [71 in which all parameters of both Gaussians are 
unconstrained, provides a better fit to the disk than the more constrained fit shown in Fig. EH 

The bottom panels in Figs. EH and [7] show the offsets of the centers of the two Gaussian 
functions relative to the disk midplane. The four cuts (A-D) are also marked on this plot. 
The third panel in Fig. [7] shows the width of the two Gaussians. 


As a cautionary note, we stress here that the two-Gaussian fit is not physically motivated 
and does not capture the true complexity of the disk. The two Gaussians do not provide 
a one-to-one match to the structures an observer would identify as the main disk and the 
warped disk (or primary and secondary disks). Nevertheless, the combination of the two 
Gaussians provides an overall good fit to the entire vertical disk profile and allows us to 
trace the general radial dependence of the vertical structure, without decomposing the disk 
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to physical components. 


The Gaussian centers (shown in bottom panels in Fig. [H] and [7} reveal a NE-SW asymme¬ 
try and trace the vertically extended disk structure most prominent at ±4”. This structure 
appears as a counter-clockwise warp in the disk that is traceable to radii <1” . The structure 
seen here has been identihed as a ’’warp” in early STIS im ages (e.g. 
based on PSF-deconvolved 2003 HST/ACS images, Gol i mowski et a 


it as an inclined secondary disk. Recently, Ahrnic ct al. (1200911 and 


Heap et al.l 120001) and, 
d2006 ) have described 


Milli ct al. (2014) have 


modeled the appearance of a two-component dust disk, considering anisotropic scattering 
and high inclination (but not exactly edge-on). 

At larger separations in the SW wing - the Gaussian center(s) are again positively 
offset above the disk midplane and this opening continues beyond the field of view of our 
images. The NE wing shows a similar, but less pronounced structure. These asymmetries 
have been recognized in early coron agrap hic studies at even larger separations and referred 
to as ’’butterfly asymmetry” (IKalas fe JewittlIl995l) . 


4.2. The Warp: Separating Asymmetric Structures 


We explore the disk warp by subtracting the contribution of the ’’main disk”, based 
on the assumption that the main disk has a perfect mirror symmetry to the main disk 
midplane, following (lLagrange et al.l 2012). Therefore, we identify four quadrants in the 
image (numbered clockwise from 1 to 4, see Fig [9]) , two of which (2nd and 4th) contain large 
flux contributions from the warped disk, while the other two (1st and 3rd) are dominated 
by emission from the main disk. By mirroring the 1st and 3rd quadrants to the main 
disk mid-plane and combining these with the original 1st and 3rd quadrant images we can 
approximate the main disk. We then subtract this combined image from the 2012 data set, 
thus producing an approximately main-disk-subtracted image (see Fig. |9|). We note that 
because the main disk is not perfectly symmetric with respect to the disk midplane this 
method is only providing an imperfect separation of the warped structure. 


This subtraction highlights the disk warp and its NE-SW asymmetry. The location and 
the asymmetry are consistent with those seen in Fig. [TJ in the NE wing the warp appears 
around 2—4”, while in the SW wing it is present from 1” to 5”; not only is the emission more 
extended in the SW, but it is much more pronounced (i.e. brighter). The mirror-subtracted 
image shows a morphology that appears to be consistent with a warp caused by an inclined 
secondary disk (see Section l4~Tj) . 

In Fig. [S] we compare the radial dependence of the disk spine - as measured by the 
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Fig. 6.— Top panel: Disk vertical profiles (blue) at Cuts A, B, C, and D fit by a sum of two 
Gaussians (sum in red, components in purple). The center of one of the Gaussians was fixed 
to the disk midplane, but all other parameters were unconstrained. Second panel: A radially 
normalized disk image shows a complex vertical structure. The locations of cuts A, B, C, 
and D are marked. Bottom panel: The vertical offset of the first- and second Gaussian fits to 
the vertical disk profiles as a function of radius. The plot reveals large NE-SW asymmetry, 
mainly introduced by the presence of the warp (within 4”). The first gaussian fits are the 
darker symbols. 
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Fig. 7.— Top panel: Same as in Fig. [61 but all parameters of the two Gaussian 
unconstrained and the radially normalized disk is shown in a 4 xvertical stretch. 
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Disk Spine Position: Observations and Models 



Distance [au] 


Fig. 8.— Comparison of the observed and modeled vertical position of the disk spine (bright¬ 


est point). Shown are composite (two-disk) models by Milli et ah ( 20141 ) and Ahrnic et ah 


( 200 9). The bine symbols show the peak position of a single vertical Gaussian function fitted 
to the STIS observations presented in this paper, an analogous measurement to the modeled 
spines shown. 
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Fig. 9.— Result of the subtraction of the quadrants without warp from the two quadrants 
with warp. The warp structure shows a clear NE-SW asymmetry. The image is rotated by 
60.9° counter-clockwise from the standard north up orientation to align the disk with the 


x-axis. 
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peak po sition _of a singl e ver tical G aussi an function - to two different recent models by 
Alimic ct ah ( 2009 1 and Milli ct al. ( 2014 1. From NE to SW the STIS single-Gaussian fit 
traces a deviation to ”below” the disk mid-plane in the figure, corresponding to the SE 
direction, between 140 and 55 au, followed by a turn to the other, NW, side of the disk . 
From NE 55 au toward the star the disk spine is located ’’above” the disk mid-plane in the 
figure, corresponding to the NW plane of the disk. The disk spine displays another break 
around SW 70 au and turns back toward the disk midplane; until 130 au the disk spine is 
approaching the disk midplane. From 130 au it displays another break and leaves the disk 
midplane. 

Both models shown as comparison in Fig. [8] use a two-component disk with anisotropic 
scattering phase functions to reproduce the disk spine as seen in projection and in scattered 
light. The Ahmic et ah ( 2009h model has bee n specifically developed to match the HST/ACS 
images presented in Goli m owski et al. (120060 . which is a data set very similar in wavelength 
to the STIS data presented her e, but d oes not extend to the inner working angle shown in our 
images. The model from Milli ct ah (1201411 has been developed to explain the VLT/NACO 
L’ observations presented in that paper; the VLT paper has a similar inner working angle to 
the STIS data presented here. 


Both papers attempt to repr oduce t he wa rp ob served in /3 Pic with a combination of 


two misaligned dust disks. The Ahmic ct ah (120091 ) study uses two axisymmetric, partly 


spatially overlapping disks, each with multi-component power-law density distributions. A 
Markov Chain Monte Carlo routin e opt imizated the disk parameters to best reproduce the 
ACS observations. The Milli ct ah (2014) study used a grid-based parameter exploration. 


We compare our observations to these models as a preliminary exploration to asses 
whether the two-component disk geometry is consistent with the observed disk spine location. 
We note that while neither of the two models matches well our new, smaller inner working 
angle data, they both offer overall similar spine geometries. While our observations do not 
directly show a secondary disk, the general spine morphology supports the interpretation 
that the inner disk is inclined with respect to the outer disk, leading to a ’’warp” as seen in 
projection. We emphasize that comprehensive scattered light modeling is required for further 
interpretation of the disk, but such an effort is beyond the scope of the current work. 


5. The Planet and the Inner Disk 

In this section we explore possible connections between the disk structure and the super- 
jupiter f3 Pic b embedded in the disk. First, we compare f3 Pic b’s orbit to the inner disk 
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structure as seen in our 2012 STIS images. As of now the orbit of /3 Pic b i s close ly monitored, 
but not yet fully known. The planet was first imaged in Nov. 2003 fiLagrange et al]i2()09 ). but 
follow-up observations could not confirm its presence until 2009 ( Lagrange et alJboiO h After 


detections, its orbit has been significantly refined (e.g. 

1 - 7 " ~ ~ 

Bonnefov et al. 

2011; 

Chauvin et al. 

2012; 

Males et al. 

2014; 

Macintosh et al. 2014: 

Bonncfoy et al. 

2014 

). Although the orbit is 


well constrained after 2009, uncertainties are very large in its position prior to 2009, based 
on a single detection. 


Here we adopt results from one of the latest orbital analyses of /3 Pic b, wh ich includes 
all published and several yet unpublished VLT/NACO observations (Lagrange et al. 2014), 
12 positions in total. The planet’s orbit has been fitted with a Markov Chain Monte Carlo- 
based (MCMC) optimization of a Keplerian orbit. With less than one quarter of the planet’s 
orbit followed accurately, multiple solutions are possible. Table [3] summarizes three different 
solutions that are consistent with the existing astrometry. The three solutions differ in the 
following ways: In Solution 1 the eccentricity was fixed to zero; Solution 2 represents the 
most likely solution from the MCMC program; Solution 3 is a high-eccentricity solution far 
from the peak of the distribution, i.e. less likely than the others, but still consistent with 
the current set of observations. During the preparation of this manuscript an independent 
orbital model ing based on exis t ing V LT/N ACO and new Gemini GPI observations has been 
presented bv Macintosh et al.l (120141) and iBonnefov et al .1 (120141) : we note that all three of 
our solutions are fully consistent with the best fit solution and its uncertainties given by 
these papers. 


In Figure [10] we plot the best-fit orbit of the planet (Solution #1) as projected on the 
plane of the sky. The plot is centered on the star’s position and we marked the planet’s 
position at a few significant epochs, including March 6, 2012, when our HST/STIS obser¬ 
vations were obtained. For the 1997 position we also plotted an uncertainty ellipse with its 
1 u semi-major axis set to match the difference in the positions predicted by Solutions #1 
and #2. 


In the left panel of Figure H21 we show the innermost section of the disk (39x39 au) 
centered on the location of (3 Pic. Black pixels mark locations covered by the coronagraphic 
wedge in every spacecraft orientation and thus without valid data. The image shown here 
is from our 2012 dataset; the earlier 1997 image has significantly larger inner working angle 
and it is less informative about the inner disk structure. 


In the right panel of Figure IT?! we compare the orbit of the planet with a contour plot 
of the inner disk surface brightness. The planet’s predicted position for the epoch of the 
observations coincides with the innermost valid pixel of our STIS data in the direction of the 
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Table 3: Orbital solutions for (3 Pic b. 


Parameter 

Solution 1 

Solution 2 

Solution 3 

a [au] 

8.6686 

9.341403 

11.18887 

P [yr] 

19.29343 

21.58237 

28.29178 

e 

0.0 

5.160x- 2 

0.163 

inc [°] 

89.269 

88.603 

88.803 

SI [°] 

-148.2210 

-148.2013 

-147.7354 

o' [°1 

169.5635 

-127.3659 

3.996902 

t p 

2002.229 

2005.643 

2013.290 

x 2 

6.26 

8.67 

5.367 
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Fig. 10.— The projected orbit of /3 Pic b (bine) and the disk midplane (red). Selected epochal 
positions of the planet at are marked with red symbols. The orbit shown here corresponds 
to Solution #1 in Table [3l The ellipse corresponds to the la positional uncertainty of the 
planet at the time of the 1997 STIS observations. 
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Fig. 11.— The warp and the planet’s orbit in (3 Pictoris. The planet’s inclination with 
respect to the disk midplane is (0.7° ±0.7, see §EJ), while the warp (vertically extended disk 
emission) is seen np to at least 4° from the midplane. 
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planet. However, no evidence for excess emission is seen in the pixel where the planet’s point 
spread function’s core should fall. This suggests that either the planet’s optical brightness 
is much weaker than that of the disk or the current best estimate for the planet’s orbit is 
slightly (~ 0"025) inaccurate and the planet may be just within our inner working angle. 


The fact that the planet is not visible inour im ages should not come as a surprise: 
Exoplanet atmospheric evolution models from Baraffe et al. (120031 ) predict a visual apparent 
brightness for (3 Pic b between V=21.2 and V=22.9 (thermal emission only, assuming a mass 


of 1 0 Mj,,p, adopting the distance of (3 Pic , for an age between 8 and 21 Mvr lOrteea et al. 


e.g. 120021; ISong et al.l e.g. 120031: lOrtega et al.l e.g. 12004; iBinks &; Jeffriesl e.g. 12014) and, based 
on simple geometric considerations, we estimate the reflected light of the planet to be fainter 
than V=26 (assuming a radius of 1.3 R j up , isotropic scattering, and a Bond albedo of 0.4, 50% 
illumination). The predicted brightness from thermal emission and reflected light is much 
lower than that of the disk at the same radius: from our 2012 image we measure 7.02 xlO -16 
ergs/s/cm 2 /A or V=16.74 for the pixel to which the planet is projected. This value must be 
a combination of scattered light from the disk, and thermal emission and reflected light from 
(3 Pic b. The fact that this pixel has very similar brightness to the surrounding pixels, which 
only include the light scattered from the disk, argues for the flux density observed at the 
location of the planet also being dominated by scattered light emission from the disk. The 
values discussed here show that the disk surface brightness overwhelms the emission from 
(3 Pic by about 5 magnitudes. Note, that the above comparison is only an approximation 
that does not correctly account for the actual bandwidth of our unfiltered STIS observations; 
nevertheless, because the choice of V-band for the comparison provides the most favorable 
planet-to-disk contrast, a more realistic estimate should End that the contribution from the 
planet is even lower. 

In Fig. ITT1 we marked the planet’s inclination (ss 0.7 ° ± 0.7° with respect to the disk 
midplane, e.g. lLagrange et al.ll2014J : Macintosh et ahll2014j) and the approximate angle under 
which the warp is observed in our STIS images (~ 4.0°). We note that the warp angle derived 
from our im ages agrees ver y well with the angle of the secondary disk structure described 
by Golimowski et ah (120061 ). The warped disk - as seen in projection - subtends an angle 
larger than the best-fit orbital inclination of /3 Pic b, suggesting that planetesimals may be 
perturbed to higher inclinations than that of the perturbing giant planet. This finding is 
consistent with the predictions of dynamical simulations of a planetesimal system influenced 


by sec ular perturbations of a planet on inclined orbits: for example, a study by (D awson et ah 


(20111) finds that the inclinations of the planetesimals will oscillate between 0 and 2 i p , where 


i p is the inclination of the planet. The fact that the warp is seen at angles 4°, would taken on 
face value, then suggests that (3 Pic b’s inclination is - within uncertainties - underestimated 
by current measurements and it may be close to ~ 2°. Alternatively, /3 Pic b’s inclination 
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may have been damped, i.e. now may be lower than in the past, which would explain a 
low present-day planet inclination with a more inclined inner disk, but would require the 
presence of another planet (Da ws on et ah 2011). However, we note that scattered light is 
an imperfect tracer of the inner disk, and the not perfectly edge-on viewing geometry and 
anisotropic scattering further complicate its interpretation. 


6. Disk Structure and Its Temporal Evolution 
6.1. Multi-wavelength View of /3 Pic 


The /3 Pic disk displays multiple radial (Sect. 13.4p and vertical (Sect. 14.ip asymmetries, 
some of which have been proposed consequences of one or two massive planets orbiting within 
the disk. In Fig. [HI] we provide a multi-wavelength-view of the disk including optical, near- 
infrared, mid-infrared, and sub-millimeter images. The top panel shows our STIS midplane- 
normalized image. The second panel from the top shows a VLT/NACO Ks-band image of 
the disk ([Lagrange et al.1 2012). For clarity we show the NACO image with two wedge-like 
sections covering the northern and southern regions directly above the star, areas that have 
been heavily contaminated by diffraction spike residuals. 

A separate paper in preparation describes the constraints from the VLT/NACO images 
and we refer to that work for a detailed comparison of the surface brightness slopes; here 
we will only focus on the disk structure and asymmetries. The Ks-band im age traces d isk 
emission in its entire field of view (~ 14” in diameter). Careful analysis bv iLagrange et ah 


(2012) showed that the disk position angle is approximately 29°.3loJo> with the precise value 
depending on details of the data reduction procedure. These authors also found a slightly 
different position angle for the SW side of the disk ( 209°. 10^:1) ; these values are in good 

agreement (within la) with the 29.1° ± 0.1° value found for our STIS data (Sect. 13.ip . 

The NACO Ks-band image also shows the warp structure seen at optical wavelengths, 
although less prominently. The disk width, as seen here, reaches its maxima at ±3.5-4''0 
separations from the star . The third pan el from the top shows an 11.7 urn mi d-infrared image 
of the disk emission from Li et ah (2012). As reported in Telesco et al, (120051 ) the MIR-clump 
- located at 2'/7 or 52 au projected separation - displays a mid-infrared color significantly 
different from the surrounding disk, indicating that grains in it differ in temperature, size 
and/or composition from the particles characteristic of the rest of the disk. 

The lower two panels show 870 /jrn continuum and 12 CO 3-2 transition images obtained 
with the ALMA sub-millimeter array (Dent et al.1 2014). Both the continuum and the CO 
image show bright peaks at the SW side of the disk. Based on earlier sub-mm images 
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Fig. 12.— Left panel: The view of the inner disk (2”x2”). Right panel: The orbit of the 
planet (blue) overlaid on a contour plot of the inner disk with the red symbol marking the 
position of the planet during the observations. Our images do not show evidence for excess 
emission from the planet, but allow probing the disk structure at the orbit of the planet. 
The image is rotated by 60.9° counter-clockwise from the standard north up orientation to 
align the disk with the x-axis. 






33 


STIS Mid-Plane Normalized Imac 



Distance ["] 


Fig. 13.— The multi-wavelength view of the disk reveals a major NE-SW asymmetry and 
a warp. While the NE-SW asymmetry is visible at all wavelengths (from optical through 
infrared to sub-mm), the warp (emission significantly above the disk midplane) is visible in 
the optical and near-infrared. The CO peak is also located above the mid plane and may or 
may not be associated with the warp. The left white line marks 2'.'7, the location of the MIR 
and CO clumps, and the right white line marks 375, the outer edge of the sub-nun dust clump. 
The blue dot on the top panel also s hows the approximate l ocation of the p lanet during the 


observ ations. Data from this paper, lLagrange et al.l (120121) . iLi et al.l (1201211 . and lDent et al. 


( 20141 )). The images have been rotated by 60.9° counter-clockwise from the standard north 
up orientation to align the disk with the x-axis. 
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Wilner et ah (2 0Ilf) argued that the mm peaks at ±3"5 (68 au) mark the locations of a 
planetesimal belt i n which collisions produce large grains. This picture is consistent with 
the model by Augere au et al.l (120010 and also with the new, higher resolution and sensitivity 


ALMA observations (Dent et all 2014) shown in Fig. [T2J 


Fig. [15] provides a visual summary of the key structural elements identified in the (3 Pic 
disk from multi-wavelength imaging. Table [4] summarizes the asymmetries observed in the 
/3 Pic disk at wavelengths ranging from blue optical to sub-millimeter. We also show a color- 
composite image (Fig. [T-Dl as an illustration allowing the comparison of the spatial location 
of the scattered light (blue), the ALMA dust continuum emission (green), and the ALMA 
velocity-integrated CO emission (red) in the disk. We group the observed disk structures in 
two major categories: apparently axisymmetrical and non-axisymmetrical structures. 


The apparently axisymmetrical , but not mirror-symmetrical disk structures are primarily 
those that contribute to the warped inner disk, seen in the optical and near-infrared scattered- 
light images. There is confident detection of extra emission above (SW) and below (NE) of 
the optical disk midplane (warp) at 2'. / 5-5 / /0 at all wavelengths. 


The apparently non-axisymmetrical structures are primarily those that reflect the major 
SW-NE asymmetry: the MIR and CO clumps present in the SW disk, the surface brightness 
differences between the SW and NE disks, and the different ALMA dust continuum levels. 
Specifically, there is evidence that the SW side of the disk is brighter than the NE at optical, 
near-infrared, mid-infrared and sub-millimeter wav elengths. New ALMA observations also 
show a similar asymmetry in CO gas (IDent et ah 2014). In addition, new Herschel/HIFI 
o bservat i ons also arg ue for a higher abundance of C II gas in the SW wing than in NE wing 
( Catald i et ah 20131 ). The MIR and mm images also argue for a different dust population 
above the disk, indicating recently released smaller grains (in the MIR) at a projected sepa- 



Fig. 14.— Color composite image of the j3 Pictoris disk. The scattered light STIS image 
is shown in blue, the 1.3mm ALMA dust continuum image is shown in green and the CO 
velocity integrated line emission is shown in red. The ALMA data are from 
(l2014h . 


Dent et al. 
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ration of 52 au and l arger grains th at trace the planetesimal belt that released them at 68 au 
flTelesco et al.l 120051: iLi et al.l 120121) . Our images also argue for a warp at this location and 
with the SW warp being brighter than the NE one. Finally, we find that the disk surface 
brightness slope in unbroken between 0"4 and 2 / .'0, indicating no significant changes in the 
disk surface density at radii adjacent and exterior to the most likely semi-major axis of the 
planet. 


6.2. Temporal Changes and Comparison to Keplerian Motion 

We now search for changes between our two images spanning 14.5 years to explore any 
temporal evolution that may be observable in the disk at optical wavelengths. We followed 
two different approaches to compare the images: In the first approach we aligned the im¬ 
ages based on the stellar position and position angles identified during the reduction process 
(based on the location of the diffraction spikes and the spacecraft pointing information). 
In the second approach we searched for the optimal offsets and 5 y) and rotation be¬ 
tween the two images by minimizing the residuals in the subtracted image. Specifically, we 
carried out a grid search covering the relative position angle range +1° to —1° (with steps 
of 0.005°) to identify the relative rotation angle that minimizes the sum of the absolute 
subtraction residuals between the two images at stellocentric separations >3.5”. For each 
relative position we cross-correlated the images to identify the optimal offsets. This method 
has identified 0.0° relative position as the optimal position with offsets of about 1 pixel. 
Tests of the second method via the introduction and recovery of artificial image rotations 
and offsets revealed that while the rotation angle is reliably identified with an accuracy of 
0.1° or better, the cross-correlation-based image offset determination’s typical uncertainty 
was ~1.5 pixels (probably set by the broad structures in the images and the noise level). 
Given that our expected positional uncertainty from the image reduction is more accurate 
than the cross-correlation based centering, we opted to adopt the first approach outlined 
above, but note that the second approach also resulted in consistent results. 

In the 14.5 year period separating our Epoch 1 and 2 images the STIS detector response 
has changed slightly. The temporal evolution of the detector’s response as a function of 
wavelength has been monitored closely the Space Telescope Science Institute since STIS’s 
activation. The evolution is described by Krist et al. (2013, STIS Instrument Science 
Report 2013-03). Standard star observations show that in the blue (filter G430L) STIS’s 
sensitivity has declined by 2.8%, while in the red (filter G750L) its sensitivity has declined 
by 2.5%. We used the Fphot,\ keywords in the FITS headers of the two datasets to convert 
our instrumental count rates into units of flux density; these conversion factors reflect the 
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Table 4: Summary of the asymmetries observed in the (5 F , ic_disk. A visual summary of 
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sensitivities calculated for the two different datasets given the unfiltered CCD throughput 
curves and assuming a flat source spectrum. Given the accuracy of the calibration and the 
low wavelength-dependence of the sensitivity change we estimate an uncertainty of ~0.3% 
in the F Photx factor. 

Figure [16] shows the Epoch 2 (2012) STIS image in a logarithmic scale (upper panel) 
and the ratio of the Epoch 1 and Epoch 2 images (lower panel). Naturally, the ratio image is 
sensitive to the signal-to-noise ratio of the images divided. The ratio image is dominated by 
noise above and below the disk midplane, and at the inner and outer edges of the images. The 
structures seen in the inner disk (<275) are dominated by PSF residuals which are different 
in the two images and result in noise. The ratio image, however, provides information on 
changes in the disk properties between 275 and ~13”. At these separations the ratio is close 
to 1.00, with variations typically below 3%. The overall morphology of the ratio image does 
not suggest that any structure within the disk midplane has changed beyond these levels 
over our 15 year baseline. 


Figure [171 shows a magnified image focusing on the 370-670 region southwest from the 
star, extracted from the larger ratio image shown in Fig. OR] This region covers most of 
the SW warp, the CO clump reported in Dent et ah (2014). The upper panel of this figure 
also shows the velocity-integrated ALMA CO emission overlaid in contours. In order to 
improve the signal-to-noise ratio we binned the image by 2x2 pixels using a flux-conserving 
1DL algorithm (frebin). The left edge of the image (~ 370) is still influenced by the PSF 
residuals; however, the rest of the image provides a higher quality measurement of the disk 
brightness evolution. This entire section of the image is within ±2% of 1.00, the level 
corresponding to no change. Given the small uncertainties in our photometric calibration we 
conclude that there is no evidence for changes in the disk surface brightness between 1997 
and 2012 at levels 3% or higher at the disk midplane or above at the location of the warp, 
the MIR clump, or the CO/sub-mm continuum clumps. 


In Fig. [18] we compare the projected Keplerian motion for circular orbits over our STIS 
baseline (14.5 years) as a function of the orbital semi-major axis (x-axis) and the ratio of the 
projected separation and the semi-major axis (y-axis). For example, the upper edge of the 
plot (y=l) shows the projected motions if the planet is seen at quadrature; smaller y values 
correspond to narrower observer-star-orbital position angles. The red-shaded contours show 
the projected motion in arc seconds, ranging from 071 to 079 in the figure. Yellow contours 
mark the x/y combinations that correspond to 37 , 470, and 570 separations from the star. 
Based on the comparison of our Epoch 1 and Epoch 2 STIS images we excluded changes 
in the disk mid plane (±0725) greater than 3% at separations 370 to 570. The resolution 
elements of the STIS image is approximately 2 pixels or 07101. This value is smaller than 
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the smallest projected Keplerian motion possible within the bounds of Fig. [I8l 

Therefore, we conclude that any point source on a Keplerian orbit in or near the disk 
mid-plane contributing more than 3% to the surface brightness distribution to the disk 
surface brightness, with projected separations between 3"0 and 5'.'0, would have led to a 
detectable change in our two-epoch comparison. This area covers part of the disk warp 
as well as the SW extension of the dust continuum seen in the ALMA images. It does 
not directly probe, however, the peak of the ALMA-detected integrated CO clump and the 
MIR-detected dust clump located at 2 , . , 7 projected separation. 


6.3. Temporal Evolution Constraints on Extended Disk Structures 


In this section we use a simple approach to translate the upper limits on change in 
disk surface brightness (see 06.21 and Fig. [RH) to upper limits on the time variations in 
the disk density structure considering a complicating factor, the disk’s edge-on orientation. 
We synthesized some images of a generic dust cloud around (3 Pic assuming an isotropic 
scattering phase function. The model images contain azimuthal structures in the shape 
of simple cosine perturbations on the dust density, expressed in wavenumbers (m). These 
structures rotate and shear with the Keplerian shear. We compared these models the region 
from 3-6” from the star on the SW side where our data is the best. 


The model has an unperturb ed dust den sity with a radial power law index of —1, and a 
vertical distribution given by the [Kelsall et ah (119981) model of the solar zodiacal cloud. To 
simulate the perturbations, we multiplied the density distribution by the function: 


1 + Ax cos(m x (A hei - Q x t) - 9 ripp i e ) (1) 

where A is the amplitude of the ripple, A^i is the stellocentric longitude, m is the 
azimuthal wavenumber, t is time, fl is the local Keplerian angular frequency, and 9 r i PP i e is 
the phase of the wave at time zero. This representation is equivalent to a cosine function 
that winds up with the Keplerian shear. 

We created images at time —7.5 and +7.5 years later (via t, see Fig. [THD and took the 
ratio of the fluxes in a slice through the midplane. We then repeatedly rotated the disk 
slightly (via 9 ripp i e ), calculated the r.m.s. of the ratios to average over rotation angles. We 
repeated this process for each of several angular wavenumbers and plotted slices to generate 
predictions for the signal as a function of radius. 

Only the m — 1 wavenumber yielded any substantial signal in the 3-6” region for 
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amplitudes A < 1. Figure H9l shows the rotation-averaged rn — 1 model compared to the 
ratio of (3 Pic images. A dashed line shows the 3% uncertainty in the ratio. The m = 1 model 
in this figure has an amplitude of 50%, i.e. A = 0.5. This 50% perturbation is marginally 
consistent with the detected change in the ratio. An rn — 1 perturbation of this kind might 
be produced via an ava lanche of dust from a single localized event - a large collision, for 
example (see Krai et ah 2014). 


Our simple model shows that our tight upper limit on the changes (< 3% at 3” and 
beyond) requires a disk that is devoid of high-amplitude, low-order angular modes and small, 
isolated structures; instead, f3 Pic appears azimuthally homogeneous at the radii probed. 


We point out that the methodology described here can be used for future multi-epoch 
observations of disks and demonstrate that sensitive multi-epoch observations of nearby 
debris disks can provide powerful constraints on the azimuthal structure of the disk. 


Our current observations are limited by the accuracy of the lower-quality first-epoch 
STIS image (1997). For the well-sampled disk regions (mid-plane, beyond 370) our com¬ 
parison yields an upper limit of only 2-3%, which is the most accurate such comparison 
yet. Interestingly, even at these large separations and long periods the 14.5-year baseline 
results in large and easily-detectable projected Keplerian motions (see Fig. [18]) . The preci¬ 
sion reached shows that further follow-up observations separated by just 3 years from our 
2012 images would yield images capable of probing orbital motions at 277, the location of 
the MIR- and CO-clumps. Similarly, image pairs separated by only 3 years, taken by the 
James Webb Space Telescope or the ALMA sub-millimeter interferometer, should be able to 
sensitively probe Keplerian motions in the disk. 


7. Possible Origins of The Disk Structure in the (3 Pic System 

We now contrast disk structure models with the high-quality multi-wavelength images 
(§H]) that cover the j3 Pic disk and with the improved orbital fits to j3 Pic b (§5]). We also use 
simple models to interpret our upper limits on temporal variations i H6.2H and demonstrate 
the future potential of such measurements. 


7.1. Planet-induced Disk Structures 


Multiple models have been proposed t o exp lain the asymmetries of the (3 Pic disk. 
Mouillet et al.l (119971) and lAueereau et al.l (120011 1 introduced a dynamical model that ac¬ 
counts for many of the observed properties of the /3 Pic disk and successfully predicted the 
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presence of a massive planetary companion in the disk. This model assumed a giant planet 
at an orbit inclined by 3° with respect to a planetesimal belt, which extends up to ~120 au. 
A collisional cascade within this planetesimal belt replenishes the dust grain population and 
the smallest grains (with sizes below the blow-out size) are efficiently moved by radiation 
pressure from the planetesimal belt to longer orbits. These small grains are responsible for 
the very large scattered-light disk seen beyond ~120 au. 


In this model the warp in the inner disk is a direct result of the gravitational perturbation 
of the planetesimal belt by the giant planet. The inclined orbit of the giant planet forces 
the precession of the planetesimals’ orbit, which leads to a warp. The outer radius of the 
warp will slowly grow in time, as the effect of the perturbations gradually accumulates for 
planetesimals on longer orbits. In addition, their simulations show the small dust grains 
generated by a collisional cascade in the warp are then subsequently blown out and will 
provide an extended, asymmetric structure, fully consistent with the prominent butterfly 
asymmetry. 


The disk structure predicted by the Auger eau et al. (120011 ) model continues to show an 
excellent agreement with the axisymmetric disk structures in f3 Pic , including the radial 
surface brightness profile (§ 13. 2J) , the radially-normalized vertical disk profile (§ 14. Tj) and 
the disk warp (§ 14.2[) , as well as the large-scale structure and the inner disk structure (§ ED- 
The model has only two assumptions: a planetesimal belt extending out to 120-150 au and 
the presence of a giant planet on an inclined orbit. The latter prediction has been verified 
throug h the detection of 0 Pic b and the former is supported by the submm observations 
(SMA: Wilner et al. 2011 . ALMA: Dent et al. 2014 1. 


7.2. A Simple Model for the Warp Morphology 

The inner disk morphology offers important insights into the disk structure and we 
employ here a simple model aiming to explain the variation in the deviations above and below 
a linear disk mid-plane with stellocentric distance characterized by the brightest points along 


Wvatt 

(2005) and its extension (see, 

Matthews et al. 

2014) 


assume an edge-on planetesimal disk with the planetesimals’ semi-major axes logarithmically 
spaced from 40-140 au, and assume that the disk surface brightness is scattered light that is 
directly proportional to the surface density of planetesimals. We placed a mass representing 
/3 Pic b on an inclined orbit (i p i = 2°) and allowed the system to evolve for 10 Myr. We 
assumed that the disk surface brightness is scattered light that is directly proportional to 
the surface density of planetesimals whose mutual collisions produce the light-scattering 



















41 


dust. This is, of course, a toy model with several important simplifications: First, the 
system parameters (semi-major axes, inclinations, etc.) have not been optimized in any 
way; second, the model only follows planetesimals and not dust grains, which are affected 
by radiation pressure. Nevertheless, we expect the model to qualitatively reproduce the 
first-order morphology of the disk and thus serves a useful way to probe the /3 Pic disk 
morphology. 

In the right panel of Fig. (2D we show the vertical location of the peak brightness at each 
radius (peak brightness profile) from this model. We vary the line-of-nodes of the planet 
(f l p i) and show the resulting peak brightness profiles compared to a peak brightness profile 
derived from the observed disk image in an identical way. We note that this representation 
is very similar, but not identical to the more detailed fitting procedure described in § 14.11 


Our simple toy model leads to three conclusions: First, we find that the orbit and mass 
of j3 Pic b and the age of the system are such that the secular perturbations from the 
planet would be expected to impose a warp on the disk at a radial location and magnitude 
that is compatible with that observed. This conclusion is in line with previous studies (e.g. 
Augereau ei ah 200 1) and confirms the trend in spite of the simplicity of our approach. 


Second, our model demonstrates the importance of constraining the shape of the warp 
in the inner region explored by our STIS observations, because that shape constrains the 
orientation of the planet’s orbital plane, specifically whether the line-of-nodes is in the plane 
of the sky (fi = 0°) or oriented along the line-of-sight (fi = 90°). If £2^=90° then the 
plane of symmetry in the inner disk would be a straight line between the warp. However, if 
fV=0° then the inner disk would be aligned with the outer disk. Early models had noted 
that this parameter affected the magnitude of the warp (iMouillet et ah 1997), but since that 
magnitude was maximized for = 90° (Fig. 6 in that paper), later models only considered 
the = 90° case. The fact that the Gaussian fits in Fig. |T] find the inner regions to be well 
aligned with the outer disk seems to favor a model with = 0° (see Fig. l2TTh However, the 
orbital fits in Table [3] favour a value closer to 90°. To explore this further would require a 
detailed model which addresses the simplifications mentioned in the first paragraph of this 
section, which is beyond the scope of the current paper. 


Third, our model also highlights the point that if the planet has an eccentric orbit then 
the same secular perturbations that result in the warp will also h ave caused the disk to 


have a tightly wound spiral structure (left panel in Fig. [2TJ e.g., IWvattll2005l and also in 


§4.3 of IMouillet et ah 1997). The spiral does not significantly affect the structure of the 
warp because the evolution of eccentricities and inclinations are decoupled for low e and I. 
However, the spiral could contribute a small brightness asymmetry to the disk. We note, 
that this spiral would be at the same radial location as the warp, since secular perturbations 
















42 


make orbital planes and pericenters precess at the same rate, which is also the same location 
as the clump. A more comprehensive modeling of the time-evolving morphology introduced 
by the projection of the spiral structure would be valuable, but beyond the scope of the 
current paper. 

With the rapidly improving orbital period and mass estimate for (5 Pic and high-quality 
multi-wavelength images at hand, a more exhaustive study of /? Pic disk dynamics is well 
motivated, but it is beyond the scope of the current paper. 


7.3. Open Questions: NE-SW Asymmetry and the Origin of a Super-jupiter 

on Inclined Orbit 


While all the axisymmetric structures observed in the disk are consistent with the pre¬ 
dictions of the models describing the dynamical interactions of a planetesimal belt and a 
giant planet on an inclined orbit, these models cannot reproduce the prominent NE-SW 
asymmetries (see Table [4] and §[6|) . This fact highlights the two fundamental open questions 
on /3 Pic : What is the origin of the NE-SW asymmetries? and What is the origin of the 
f3 Pic b super-jupiter’s inclined orbit? While our observations do not directly help to answer 
the latter question, they provide some constraints on the different scenarios proposed to 
explain the former. 


The mid-infrared observations by Te lesc o et al. (1200511 and Li et al. (2012) show that 
the dust grain population in the bright dust clump in the SW wing is different (possibly 
smaller grains) than those typical to the rest of the disk. These authors argued for the 
possibility of a single, recent collision injecting copious amounts of dust in the system. The 
recent ALMA images by Dent et al. (2014) show that the SW clump not only has a different 
dust population, but also contains gas-phase CO. Given that the UV-photodissociation 
li fetime of gas-phase CO by the interstellar UV field in an optically thin disk is ~120 years 


Visser et al.l (120091) . the observed CO gas mu st h ave been released recently, consistent with 


a recent major collision (ITelesco et al.l 120051 and iDent et al.l I2014J) . The estimated mass in 
the SW dust clump detected by ALMA argues for a Mars-s ized parent body, assuming that 
about 10% of its mass has been released as d ebris (IDent et al.ll2014j) . A similar recent major 
collision has been proposed by S tark et al . (2014) as one of two possibilities to explain the 
asymmetries observed in the HD 181327 debris disk. To avoid the improbab ility of having 
witnessed such a massive collision only 120 years ago in /3 Pic , Jackson et al. (2014) pointed 
out that the asymmetric dust distribution resulting from such an event can last up to 0.5 Myr 
if the dust clumps mark the location where the collision occurred, since debris released by 
the collision will have a variety of orbits but the orbits will converge at the collision point. 
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However, this scenario would predict that the dust clump’s location is constant (i.e. not 
on Keplerian orbit), which is marginally inconsistent with the tentative detection of orbital 
motion of the clump Li et al. (2012). 


An alternative scenario has been proposed for ft Pic based on work by IWvattl (120031 ) 
and W yatt (120061 ). In this picture a planet - during its outward migration - captures plan- 
etesimals in mean motion resonances. The frequent mutual collisions of these planetesimals 
give rise to dust clumps that trace the planetesimal population and thus produce resonant 
structures that co-rotate with the planet. This scenario predicts a >10 M Earth planet at the 
inner edge of the dust belt (with half the orbital period, as that of the be lt). This scenario 


has been proposed for ft Pic both by iTelesco et al.l (120051) and iDent et al.l (120141) 


Perhaps the most directly testable key difference between these two scenarios is the SW 


dust clump’s orbital motion or its absence: if the clump is on orbit (as suggested by Li et al. 


2012) then this would argue for it being a resonant structure. If its location is constant, the 
clump would likely be a result of a recent collision. 

It is important to point out that neither of the above scenarios for the origin of the SW 
emission invoke the dynamics proposed to explain the warp. Yet, the warp is at exactly the 
same radial location as the clump, which may be either a coincidence or there is a causal 
connection. There may be scenarios in which the warp aids the formation of the clump. 
For example, if ft Pic b is on an eccentric orbit then there will be a spiral pattern at 
the location of the warp. The resulting brightness asymmetry is much smaller than that 
observed. However, we also know that this is al so t he l ocati on where orbits of debris that 
were previously not overlapping start to cross (IMustill fe Wvattl 2 0091 ). which might increase 
the chances of seeing a recent collision at that location. For the second scenario there is not 
a causal connection, but the spatial coincidence is not surprising, since the orbit of the outer 
planet (that with the resonantly trapped planetesimals) would be aligned with ft Pic b and 
would also have aligned the resonant planetesimals with its orbit. 


8. Summary 

We present new and re-reduced archival HST/STIS coronagraphic optical imaging of the 
ft Pictoris debris disk. Our data provide the yet highest signal-to-noise and smallest inner 
working angle optical images of the disk. Based on our images we characterize in detail the 
disk structure and compare the optical images to published images taken at wavelengths 
ranging from the near-infrared to mm. 

The key findings of the paper are the following: 
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1) Our high-contrast images provide a continuous coverage of the inner disk covering 
radii from 0"35 to 13” in unfiltered light. The inner working angle is about 2x smaller than 
that of previous optical images and allows us to image the disk at a location and time where 
the j3 Pic b planet was present in 2012. 

2) We identify radial regions with constant radial surface brightness slopes. We show 
that corresponding regions in the northeastern and southwestern wings have similar, but 
slightly different surface brightness slopes, indicating an asymmetry in the disk. 


3) We find that the surface brightness slope between 0"4 - 2 / /0 is constant, arguing 
against any significant changes in the disk structure at or adjacent to the most likely semi¬ 
major axis of (3 Pic b. 


4) We present the first optical images of the inner structure of the main disk and its 
vertical extension. While our im ages d o not show the presence of a separate secondary disk as 
suggested by Golim owski et al. (120061 1. preliminary comparison to two-component scattered 
light disk models suggest that our observations are consistent with the warp caused by the 
projection of an inner, inclined disk onto the outer disk. 


5) NE-SW asymmetry: we confirm that the two wings of the disk have different radial 
surface brightness slopes and that the SW wing is brighter in the inner 8"0 but fainter beyond 
that than the NE wing. We show that the same asymmetry extends to the innermost disk, 
down to at least (y/5. 


6) The angle of the disk warp (seen in projection) is approx. 4°, significantly larger than 
the best-fit inclination of f3 Pic b’s orbit (0.7°). 

7) Careful comparison of ST1S images obtained over 15 year baseline (in 1997 and in 
2012) shows no difference greater than 3% in the disk surface brightness at the locations of 
the disk warp part of the CO and sub-mm continuum clumps (SW 3-5”). Similarly, we do 
not detect any difference between the radial surface brightness profiles. 

8) We compile all disk asymmetries seen in wavelengths ranging from optical through 
mid-infrared to sub-nun wavelengths. We divide the asymmetries in two groups: apparently 



collision or the presence of plantesimals on resonant orbits with a yet unseen planet at 
~80 au. 


9) We show that over a few year baseline projected Keplerian motions are sufficiently 
large to allow detection with existing and near-future facilities. We argue that, in particular, 
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HST, JWST, and ALMA multi-epoch observations will provide powerful constraints on the 
azimuthal structure, dynamics, and disk-planet interactions in the ft Pic and other nearby 
debris disks. 

Our work leads to two fundamental questions required to explain the properties and 
origin of the f3 Pic system: How did super-jupiter (3 Pic b ended up on an inclined orbit? 
and What is the origin of the NE-SW asymmetry? We discuss different scenarios and show 
that time-resolved studies of Keplerian motions in the ft Pic can be a powerful way to 
discriminate these. 
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Charles Proffitt for their dedicated support of this program. We thank Jean-Charles Augereau, 
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the Association of Universities for Research in Astronomy, Incorporated, under NASA con¬ 
tract NAS5-26555. MCW is grateful for support from the EU through ERC grant number 
279973.We acknowledge support from the French National Research Agency (ANR) through 
the grant ANR10-BLANC0504-01. 
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Sketch of the Beta Pictoris System 
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Fig. 15.— Key structures in the j3 Pic system, as derived from multi-wavelength imaging. 
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Log-scaled Epoch 2 
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Fig. 16.— Upper panel: Logarithmically scaled Epoch 2 disk image. Lower panel: The 
ratio of Epoch 1 to Epoch 2 image probes disk surface brightness evolution between 1997 
and 2012. The ratio image inside 3"0 is dominated by PSF residuals from the Epoch 1 
image. Between 3"0 and 5'.'5 in the disk midplane no significant changes are observed in the 
disk structure. The images have been rotated by 60.9° counter-clockwise from the standard 
north-up orientation to align the disk with the x-axis. 
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Log-scaled Epoch 2 and ALMA Integrated CO Contours 
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Fig. 17.— Top: STIS direct image of the disk in logarithmic scale and velocity-integrated 
CO intensity from the ALMA observations by Dent et al. (2 0141 ). Bottom: Relative change 
between Epoch 1 and 2 in the SW wing. In the disk mid plane between 3-6” our mea¬ 
surements do not show changes greater than 3%. The images have been rotated by 60.9° 
counter-clockwise from the standard north-up orientation to align the disk with the x-axis. 
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Projected Keplerian Motion in Arcsec between Epochs 1 and 2 
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Fig. 18.— The red-shaded contours show the difference in the projected location of a particle 
(in arc seconds) at a Keplerian orbit around f3 Pic between our Epoch 1 and 2 observations. 
The yellow curves show the curves corresponding to 3”, 4”, and 5” projected separations 
from (5 Pic . The STIS pixel scale is (y/0507 and the predicted projected orbital motion is 
greater than a resolution element over the entire plot. Our observations place a 3% upper 
limit at the difference in the mid-plane disk surface brightness over 15 years, effectively 
excluding the presence of single dust clumps contributing at this level to disk above the 3” 
yellow curve. 
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Fig. 19.— Face-on images of the m=3 mode, at time 0 and 15 years later, in logarithmic 
scaling. For the plot of the predicted temporal evolution signal shown in Fig. [20] edge-on 
images were used. 



ratio 


-51 - 

1.06 

1.04 

1.02 

1.00 

0.98 

0.96 

3 4 5 6 

arcsec 

Fig. 20.— Change in surface brightness between Epoch 1 and 2 for the angular wave number 
(m = 1, assuming 50% perturbations, i.e. 4=0.5). Our upper limit of 3% beyond 3” argues 
for absent or very low-amplitude low-order wavenumbers. 
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Fig. 21.— Left: Our simple model for disk-planet interactions predict a tightly wound 
spiral in the density distribution of the disk if the planet is on an eccentric orbit. Right: The 
vertical location of the peak disk brightness at each radius as a function of the orientation 
of the linc-of-nodes (in color) and the same observed quantity (thick solid line). 
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